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Abstract. Quantum spin-1 chains may develop massless phases in presence of Ising-like and single-ion 
anisotropies. We have studied c = 1 critical phases by means of both analytical techniques, including a 
mapping of the lattice Hamiltonian onto an 0(2) NLctM, and a multi-target DMRG algorithm which allows 
for accurate calculation of excited states. We find excellent quantitative agreement with the theoretical 
predictions and conclude that a pure Gaussian model, without any orbifold construction, describes correctly 
the low-energy physics of these critical phases. This combined analysis indicates that the multicritical point 
at large single-ion anisotropy does not belong to the same universality class as the Takhtajan-Babujian 
Hamiltonian as claimed in the past. A link between string-order correlation functions and twisting vertex 
operators, along the c = 1 line that ends at this point, is also suggested. 
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1 Introduction 

One-dimensional quantum spin systems have been exten- 
sively studied, since when [1], twenty years ago, Haldane 
argued that half-integer spin Heisenberg antiferromagnetic 
(AF) chains have no spin gap and are quantum criti- 
cal, while integer spin chains are in the so-called Haldane 
gapped phase and exhibit correlation functions that decay 
exponentially with a finite correlation length. This sce- 
nario has then been confirmed both numerically [2] and 
experimentally [3]. 

We will consider here instead the following spin-1 Hamil- 
tonian 

L 

H — {^S^S^j^i + SjS^^i + XSjSj_^_l + £'(S'|)^} , 

(1) 

which includes both an Ising-like and a single-ion anisotropy 
term, with coefficients A and D respectively. In our cal- 
culations we impose periodic boundary conditions (PBC): 
Sl+1 = Si. The inclusion of such terms is relevant for a 
better understanding of experimental compounds such as 
NENP and CsNiCla ([4] and references therein). 

Even if no exact solution of (1) is available (except for 
spin-1/2 [5]), it is known from numerical studies [6] that 
the inclusion of anysotropy terms can drive the system 
away from the Haldane phase toward other phases, some 
of which are critical. A first theoretical study which de- 
scribes the critical properties of the model can be found 
in [7] where, in general a spin-S* chain is mapped onto a 



system of 25" coupled spin-1/2 chains and analyzed with 
bosonization techniques. More recently, a thorough ana- 
lytical study of this model has been presented in the sem- 
inal work [8], where the emphasis is put on the physical 
properties of the massive phases. 

The basic structure of the ground-state (GS) phase di- 
agram of (1) appears to be well understood [9] and shows 
a rich variety of phases. For high values of D the system 
is in the large-D phase (D) consisting of a unique GS with 
total magnetization S^^^ = separated by a gap from the 
first excited states which lie in the sectors S^^^ = ±1. For 
large positive values of A we have a twofold degenerate 
AF Ising-like phase (I). For not too large D and A, these 
two phases are separated by the Haldane phase (H), which 
includes the isotropic 0(3)-symmetric point. It is charac- 
terized by non vanishing string-order parameters 



k-l 



lim 



5;exp iTT J2 Sn\S"k 



(2) 



n=j + l 



with a — {x, y, z}, first introduced by den Nijs and Rom- 
mclsc [10]. The first excited states are magnons at the 
boundary of the Brillouin zone (BZ) carrying total spin 1. 
These three massive phases can be distinguished on the 
basis of a hidden Z2 x Z2 symmetry [8], which is fully 
(in the Haldane phase) or partially (in the AF Ising-like 
phase) broken and whose order parameters are given by 
Eq. (2). Also, the H-D and the H-I transition lines meet 
at a tricritical point, for I? ~ A ~ 3. For larger values 
of the parameters, the Haldane phase disappears and the 
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Hne D \ represents a first order transition between the 
large- and the Ising-like phases. The remaining portion 
of the phase diagram, for A < 0, consists of a ferromag- 
netic GS for large |A| and of two critical gaplcss phases, 
XYl and XY2, different for having, for finite-size systems, 
first excitations carrying S^^^. = ±1 and S^^^ = ±2 respec- 
tively. 

In this paper, we will investigate the critical properties 
of the Hamiltonian (1), and in particular we will concen- 
trate on the Haldanc/large-i? (H-D) transition line and on 
the XY2 massless phase, which turn out to be described by 
c = 1 Conformal Field Theories (CFT). Wc will tackle the 
problem both via analytical tecniques, by finding the CFT 
Lagrangian that describes the model in the low-energy 
continuum limit, and via a multi-target Density Matrix 
Renormalization Group (DMRG) algorithm, which allows 
for accurate calculation of excited states, to be compared 
with the operator content of the CFT. We will try also to 
clariiy some controversial aspects discussed in the litera- 
ture. 



2 The 0(2) NLaM on the H-D transition line 

The H-D transition line has been located numerically us- 
ing the twisted boundary method in [9] , where it has also 
been pointed out that this represents a second order phase 
transition described by a c = 1 CFT, in accordance with 
[7]. Here we will describe a mapping of the lattice model 
(1) along the H-D transition line onto an 0(2) nonlinear cr- 
model (NLcrM) , that establishes a connection between the 
coupling constants D, A of the discrete model and those 
of the continuum Gaussian theory, namely the spin-wave 
velocity v and the compactification radius. This will allow 
us to make quantitative predictions. 

The partition function for Eq. (1) in a path-integral 
representation which makes use of spin coherent states is 
given by [11]: 

Z = y[i)r?]exp|i5^a;[r?(j,r)]-^''drF(r)| , 

(3) 

where the vector operator Sa{j) has been replaced by the 

classical variable sf2a{j,T) and ui[ila{j,T)] is the Berry 
phase factor. In a semiclassical approach, we can expand 
sf2a{j,T) about the classical solution which, for D > A — 
1, is a planar state where the the unit vectors ilaii-.T) 
are Neel ordered in the a;y-plane: fl{j,T) = (cos(6'o + 
j7r),sm{0o+jn),0). Hence we make the Haldane-like ansatz: 



taking into account that the Berry phase is given by 



f2{j,T) = {-iyn{j,r)\ll- ^2 



+ z 



(4) 



where n{j,T) = e'^^-'''^) e 0(2)^^, z is the unitary vector 
(0,0,1), and the fluctuation field is supposed to be 
small. Expanding H{t) up to quadratic terms in and 



is ^ uMj, T)]=iJ2 [ dT I (j, r) a 



0{3,t), (5) 



we can now integrate out the fluctuating field. If we treat 
then 6{j, r) as a slow-varying variable, in the continuum 
limit we end up with an effective 0(2) NLcrM in the field 
6 that, after the rescaling = 6/ sjg, can be rewritten in 
the standard form: 



-0(2) 



where 



5= -v/2(l + -D + A); t; = 5^2 (! + £> + A) 



(6) 



(7) 



This is a free Gaussian model [12], with a bosonic field Q 
compactified along a circle of radius 1 / y^, which describes 
a CFT with central charge c = 1 and with primary fields 
(vertex operators) Vmn of scaling dimensions given by 



m,n £ 



(8) 



where K = n/g and n and m are, respectively, the winding 
numbers of and its dual field, The latter turns out 
to be compactified on a radius y/g/2TT. Moreover, it is not 
diflacult to see that m = Sf^j, and hence it is a conserved 
quantity of the model. Indeed, from the NLcrM approach, 
it follows that S^'(x) l(x) — drO/vg so that 



Hot 



/dx = [ dx = ^— 
^^v^ J V9 27r 



m 



(9) 



We recall that the values K 



,1,2 correspond to the so- 



called [13] self-dual (SD), free Dirac (FD) and Berezinski- 
Kosterlitz-Thosuless (BKT) points respectively. The scal- 
ing dimensions (8) fix also the (non universal) critical 
exponents of the correlation functions. For instance the 
transverse spin-spin correlator decays according to: 



{S+{0)S-{x)) « (e'^(o^e-'^W) qc 



(10) 



where r] = 2rfio = g/2'jT. 



Let us turn out to numerical results on the Hamilto- 
nian (1). We used a DMRG algorithm [14] for finite sys- 
tems, which had to be customized for the convergence of 
not only the GS but also of several excited state energies. 
More specifically, using the so-called thick-restart Lanczos 
algorithm [15] we could target up to eight states in a given 
sector of Sfgf . Then, the density matrix is built by averag- 
ing (with equal weights) the density matrices associated 
with these states. Shortly, we will call this a multi-target 
DMRG. If necessary, the correlation functions on the GS 
can be computed at the end of the finite-size iterations 
(three in our cases). Since the DMRG is known to be a 
good tool to investigate systems with a relatively short 
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correlation length, in order to study critical phases we had 
to handle the DMRG data with finite-size scaling (FSS) 
techniques. For the methodological aspects of our DMRG 
procedure, and fitting of data, we refer to a forthcoming 
paper [16]. 

First of all, it has been necessary to locate the H-D 

transition line with great precision. Initially, we have fixed 
some representative values of A and let D vary across the 
phase boundary by small increments. Then, we have re- 
fined the location of the critical points Dc{X) according to 
FSS theory [17], using M = 400 DMRG states for chain 
lengths L ranging from 10 to 50. The so obtained values 
are very close to the ones calculated in Ref. [9] using the 
twisted boundary conditions method and exact diagonal- 
ization limited to X = 16 sites. 

It is well known [18] that in a CFT of a finite size 
system of length L. the GS energy density depends on the 
central charge and converges to its thermodynamic limit 
as 



The excited state energies are instead related to the scal- 
ing dimensions (8) by 

Emn - Eoo = ^ {d,nn + r + f) , (12) 

where r and f are positive integers that label the sec- 
ondary states of a Verma module [19]. 

In a numerical approach, Eq. (11) is the starting point 
to identify the correct CFT for a given critical point: accu- 
rate calculations of Eqo at various L give in fact a best-fit 
of Coo and of the product cv. Then one has to select a 
number of excited states that become critical, i.e. degen- 
erate with the GS, in the limit L oo. In our problem 
the energies of the lowest states have been calculated for 
different values of which is the only quantum num- 
ber that can be fixed within the DMRG algorithm. If the 
hypothesis of an underlining c = 1 Gaussian CFT is cor- 
rect, all these energies should approach zero as straight 
lines as functions of 1/L, according to Eq. (12) with dmn 
given in Eq. (8). For the critical point (A = 0.5, D = 0.65) 
this is shown in Fig.l, where we have also indicated the 
(quasi-)degeneracy of each state, as observed numerically 
for chains of any length. The numerical data (points) of 
Fig.l have been best-fitted with lines, whose slopes rf^"""*) 
are reported in the second column of Table 1. 

In order to compare this spectrum with the theoretical 
predictions, it is necessary to find the values of ?;(""™) and 
j^inum)^ that identify the Gaussian model from which the 
observed levels originate. Eq. (8) shows that, \i K > 1, 
the first excited state in the m = sector should be 
a doublet, corresponding to rn = n = and (r, f) = 
(1,0), (0, 1). Numerically, we do indeed find a doublet as 
first excited states in the S'f^j = sector, so that, us- 
ing Eq. (12) together with the data that come from the 
fit on the GS energy of Eq. (11), we are able to obtain 
yinum) ^ 2. 197 ±0.004 and c^""") = 1.008 ±0.003, which 
confirms the c = 1 theory behavior. Also, again from Eq. 
(8), it follows that the primary state (m = l,n = 0) 




1/L 

Fig. 1. Energy differences, divided by 2-w, plotted vs 1/L at 
the H-D transition point (A — 0.5, D = 0.65). The legend on 
the right indicates the Stot quantum number and the order 

of the targeted level within that sector, respectively. Points: 
multi-target DMRG data (with M = 405 states). 



should be the lowest one in the S^^^ = 1 sector, with di- 
mension dio = (4ir("«'"))-i. Our DMRG data yield then 
j^inum) ^ i^gQ ^ 0.004 (1). For comparison, if we plug 
the coordinates of the critical point (A = 0.5, D = 0.65) in 
the formulae (7) for v and g of the 0(2) NLctM we obtain 
^(th) ^ g(th) ^ 2.07 and if^*'') = 1.52, which confirm the 
validity of our theoretical approach. 

We can now use the so calculated and 
to obtain the scaling dimensions d!^'^^'^^ of the low-lying 
levels as predicted by (12). These are listed in the first 
column of Table 1 and are to be compared with the corre- 
sponding numerical observations given in the second col- 
umn. One can see that the differences lie within a 2% 
percent . As a final check of the validity of the NLcrM ap- 
proach, we have computed directly the transverse spin- 
spin correlation funcion, finding [16] that it decays alge- 
braically with a critical exponent r/""™) = 0.312 ± 0.002, 
in very good agreement with the theoretical value t] = 
2dw = 0.316, obtained from Eq. (10). 

We would like to stress here that our numerical anal- 
ysis shows that the spectrum of the lattice Hamiltonian 
(1) is completely exhausted by the levels of formula (12). 
This is at variance with some claims that have appeared in 
the literature in the past [10] according to which the H-D 
transition line should be in the the same universality class 
as the Ashkin- Teller (AT) model. Indeed it is known [20] 
that the critical properties of the latter are described by 



^ Throughout this paper the reported rmrricrical errors orig- 
inate from the best fits, with the exception of the velocities for 
which the indicated errors represent the spread of the quadratic 
extrapolations in 1/L of the secondaries of the GS. 
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(0,0) 
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(0.0) 
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(o.-i) 


((),{)) 
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1.589 ± 0.007 


(±2,0) 


(1,0) 




1.589 ± 0.007 


(±2,0) 


(0,1) 


1.738 ±0.004 [x4] 


1.693 ± 0.008 


(±1,1) 


(0,0) 




1.693 ± 0.008 


(±1,-1) 


(0,0) 



Table 1. Spectrum of scaling dimensions at the point (A = 
0.5, D = 0.65) on the H-D line, obtained from the scaling plots 
in Fig. 1. Here (H-D line) m gives directly the eigenvalue of 
(t Setting S'^^'^^ = d*""'"' fixes the value of JfC"^"") = 
1.580). 
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1.328 ± 0.004 [x2] 
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(0.1) 
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(0-1) 
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1.516 ±0.005 [x4] 


1.468 ± 0.007 


(±1,1) 


(0,0) 




1.468 ± 0.007 


(±1,-1) 


(0,0) 


1.695 ± 0.005 [x2] 


1.679 ± 0.005 


(±3,0) 


(0,0) 


1.750 ± 0.02 [x4] 
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Table 2. Spectrum of scaling dimensions at the point (A = 
1,D = 0.99), obtained from the scaling plots in Fig. 2. Recall 
that m = S!„t. (t Setting S'^^'^^ = d^"""") fixes the value of 



c = 1 orbifold CFT models [13] and hence should exhibit 

if-independent scaling dimensions as well. 

This result can be further checked by moving along 
the critical line, varying K. Prom a theoretical point of 
view this would correspond to acting by means of the 
marginal operator (8^6)'^. Consider for example the point 
{X = 1, D = 0.99), i.e. the D-induced transition point for 
the isotropic Heisenberg model. Numerically we have es- 
timated = 2.588 ± 0.006, c(""'"^ = 0.997 ± 0.003 
and = 1.328 ± 0.004. Again, we have a good con- 
firmation of the 0(2) NLcrM predictions, corresponding 
to = 2.45 and K^*''^ = 1.285. For this case, scaling 
plots are displayed in Fig. 2 and the scaling dimensions are 
listed in Table 2. Also, the calculated value of the critical 
exponent ^ 0.374 ± 0.003 matches the theoretical 
value r] = 2dio = 0.377. 



c = 1 CFT for spin-1 chains 




1/L 

Fig. 2. Energy differences, divided by 27r, plotted vs 1/L at 

the H-D transition point (A = 1,D = 0.99). The legend on 
the right indicates the Stot quantum number and the order 
of the targeted level within that sector, respectively. Points: 
multi-target DMRG data (with M = 405 states). 

3 Towards the tricritical point 

As just recalled, some time ago it has been conjectured 
[10] that the H-D transition line is described by c = 1 
orbifold models. It has also been argued [7] that the tri- 
critical point at which the c = 1 H-D transition line meets 
the c = 1/2 H-I transition line is in the same universality 
class as the integrable Takhtajan-Babujian (TB) Hamil- 
tonian [21], whose low- lying excitations are described by 
a SU(2)2 Wess-Zumino-Novikov-Witten (WZNW) model, 
with central charge c ~ 3/2. Since this model is equivalent 
to three free relativistic Majorana fermions [22], one can 
conclude that the two critical lines must meet at the exact 
point where the c = 1 operator content is described by a 
bosonic theory with K ~ 1 [12]. If so, moving towards the 
tricritical point on the H-D line one should find decreas- 
ing values of approaching 1 from above. Previous 
studies on the X-D phase diagram [9] give £> ~ A ~ 3 as 
a good estimate for the; location of this point. 

In the previous Section, we have seen that our NLctM 
predictions seem to be quite reliable, at least fov D < 1. 
However, the assumptions under which we have derived 
the NLcrM should be more and more valid as D (and A) 
increases, since the true GS of the Hamiltonian (1) is bet- 
ter and better described by the planar configurations of 
Eq. (4). Tweaking the values of the microscopic parame- 
ters, we find that, in our approach, the condition K = 1 is 
fulfilled for A ~ 2. Thus our theoretical scenario predicts 
that the portion of the H-D line for A > 2 must correspond 
to ii' < 1. As we will see, our DMRG data confirm this 
hypothesis and at the same time give an estimate for K 
at the tricritical point which is very close to the SD value 
1/2. 
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To check this numerically, we have considered the crit- 
ical point (A = 2.59, -D = 2.30). For such a point the 
analysis of the spectrum is more complicated than before 
because of two reasons. First, the H-D transition line is 
now very close to the H-I line and we start seeing a super- 
position of the two conformal spectra. However, the chosen 
point not being exactly on the H-I transition line, from our 
numerical data we are still able to trace and separate the 
levels of the c = 1/2 theory. The price to be paid is that we 
had to target up to eight excited states within the Sf^^ = 
sector. Second, since we are moving towards the SD point, 
wc have to take into consideration that the operators Vo±2, 
with scaling dimension dQ±2 — AK , are becoming less and 
less irrelevant (truly marginal for K = 1/2). Thus, in a 
finite size system, we might expect [23] such operator to 
induce some renormalization effects on the energies of the 
model, which might be very large as wc approach the SD 
point, where they are logarithmic. Another subtle point 
is the choice of the states from which the velocity is com- 
puted. When if < 1 we have to remember that the doublet 
with m = and n = ±1 comes first than the doublet of 
secondaries with (r, f) = (1, 0), (0, 1). Even if we have no 
direct control on how the aforementioned finite-size cor- 
rections act on this doublet, a quadratic extrapolation in 
l/L yields = 3.70 ±0.04, that seems to be quite re- 

liable in that it gives c^""™) = 0.99 ±0.01. The theoretical 
value predicted hj our NLtrM mapping is t;^*''^ = 3.43. On 
the other hand, in the estimate of X*^"'""). the problem 
can be circumvented by using the dimensions do±i instead 
of d\o as done before. In fact, the doublet with to = and 
n = ±1, which should be degenerate according to (8), is 
found to be split. Following [24], we learn that the renor- 
malization effects on the energies of these two states are 
equal but with opposite signs, so that the correct value of 
the energy can be obtained by considering their semisum. 
So doing, we get = 0.85 ± 0.01, confirming that 

we have gone closer to the SD point beyond the FD one. 

Finally, we have considered the point (A = 3.20, D = 
2.90), which according to [9] corresponds to the location 
of the tricritical point. Similarly with what we have ar- 
gued above, we find that the numerical spectrum can be 
interpreted as a superposition of a c = 1/2 and a c = 1 
CFT's. The calculation of c is complicated by the fact that 
the Ising transition is even closer, if not coincident. Using 
the same method of the previous point, our best values 
in this case are c^""'") = 1.133 ± 0.006 and = 
0.526 ± 0.007, having estimated = 4.445 ± 0.005 

against f^*''^ = 3.77 from the expression of the NLctM. In- 
terestingly enough, if we try to extract the Ising velocity 
from the secondaries of the primary state with d = 1/8 
we find the pretty close value v^^™^ = 4.35 ± 0.06. 

In Table 3 wc summarize the values of v and c for all 
the four H-D critical points discussed above and we list 
also the GS energy per site Coo and the final estimate of 
K. As anticipated, moving to the right on the H-D line 
the value of K keeps on decreasing towards the SD point 
where we speculate that this line meets the H-I one and a 
first order transition line begins. 





V 


c 


Coo 


K 


V 


(0.5,0.65) 


2.197 


1.008 


-0.908765(9) 


1.580 


2.38 


(1.0,0.99) 


2.588 


0.997 


-0.859152(2) 


1.328 


1.49 


(2.59,2.30) 


3.70 


0.99 


-0.675099(5) 


0.85 


0.870 


(3.20,2.90) 


4.445 


1.133 


-0.59132(2) 


0.526 


0.678 



Table 3. Velocity, central charge and GS energy density (er- 
rors on the last figure in parenthesis) for the critical points 
discussed in the text. The numbers arc the outcome of DMRG 
calculations with L = 16, 20, 24, 32, 48, 64 and M = 405 for the 
first two Unes and L = 16, 20, 24, 28, 32, 36, 40 and M = 400 for 
the other cases. The last two columns contain, respectively, the 
estimate of the K parameter of the effective c = 1 Gaussian 
theory, according to the numerical procedure described in the 
text and the critical exponent i/ = 1/(2 — K). 



4 Away from the H-D critical line 

So far, we have considered what happens when moving 
along the c = 1 critical line in the direction of increasing 
A. When we instead move in the opposite direction, the 
velocity v gets smaller and the compactification radius 
of grows. In particular, when the H-D transition line 
meets the A = axis, for D = 0.4 [25], we reach a BKT 
transition, for which K = 2. For A < a two-dimensional 
critical region, corresponding to the XYl phase, opens up. 

To understand what happens in the XYl region as well 
as when wc move away from the H-D line with A > 0, 
we need to consider all relevant operators, allowed by 
the symmetries, that can be generated by renormaliza- 
tion of the lattice Hamiltonian. In a Gaussian theory with 
1/2 < K < 2, formula (8) shows that the most rele- 
vant primary field not forbidden by the conservation of 
the total magnetization along the z-axis corresponds to 
(to = 0, n = ±1). It has scaling dimensions do±i = K and 
it is given by the vertex operator cos {\/4ttK<P) . 

After a dual transformation on Eq. (6) we find that 
our model has to be described by the Lagrangian of the 
sine-Gordon model 



^SG 



cos (\/4^#) . (13) 



Here, we assume that the coefficient goes to zero along 
the H-D transition line. The precise determination of (jl as 

a function of A and D would require an exact renormal- 
ization procedure [7,26], which goes beyond the scope of 
this paper. The relevant cosine term is thus responsible of 
the opening of a gap as soon as we move away from the H- 
D transition line. In passing, we point out that, once the 
value of fio±i is known, the relation v = 1/(2 — do±i) yields 
a much better estimate, as compared to the /3-function 
method [16], of the exponent that controls the opening 
of the energy gap proportional to |_D — Dc\'^ . The values 
computed with the former method are given in the last 
column of Table 3. 

The gap-generating term becomes marginal exactly at 
the BKT point if = 2, (A = 0, D ~ 0.4) [25], where the 
H-D line "fans-out" [26] into the bidimcnsional critical re- 
gion XYl. We argue that in this critical phase the effective 
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theory is the same as in Eq. (13), the only difference be- 
ing the irrelevance of the operator cos (V 'iirK^), with K 
ranging from 2, at the BKT boundary lines, towards oo 
at the ferromagnetic transition. In order to support this 
picture, we note that the NLaM approach predicts a crit- 
ical stripe enclosed between £> = — A — I + tt^/S (where 
= 2) and Z? = -A - 1 (whore X^*'*) = oo), that 
overlaps with the real XYl region in a wide portion of 
the diagram. As a numerical test, we have considered the 
point (A = —0.5, D = 0) where wc expect K^*'^^ — tt and 
y(.*-h) _ ^ Using the same DMRG procedure as above (L 
up to 40 with M = 400), we find that the low- lying spec- 
trum is again described by a purely Gaussian CFT with 
^{num) ^ 1.11535 ±0.00005, c(""'"^ = 0.9997± 0.0001 and 
j^inum.) = 3.086 ± 0.002. 

To finish our discussion on the H-D transition, we ob- 
serve that the operator content of the microscopic Hamil- 
tonian (1) docs not contain the so-called twisting operator 
^o±i/2 = cos {^ttK(1>) with scaling dimension (io±i/2 = 
K/4:. In general, this kind of operators is permitted in 
CFT [24.27], in as much as they yield half- integer con- 
formal spin, and hence well-defined correlation functions 
[12] decaying as power laws. Specifically, for Vo±i/2 the 
exponent has to be t]s,z = 2rfo±i/2 = K/2. 

This relevant operator is forbidden here only because 
of the PBC's wc have chosen and not from symmetry con- 
siderations. Nonetheless, it is interesting to note that we 
do find certain correlation functions that, both on analyti- 
cal and numerical grounds, decay with the exponent writ- 
ten above. Namely, we considered the longitudinal string 
correlator [a = z in Eq. (2)] and extracted its decaying 
exponent at the critical points of Table 3 by means of a 
proper ESS analysis [16,17] on the DMRG data with L = 
32, 48, 64, 80, 100 and M = 300. We get 4""™^ = 0.804 ± 
0.003 at (A = 0.5, D = 0.65) and 4™™'' = 0.741 ±0.002 at 
(A = 1,D = 0.99), to be compared with Vs^f'^^ = ^•'^90 
and Vs^'f^'' — 0-664, respectively. In addition, this prelim- 
inar identification is enforced by the continuum version of 
the string correlations in the framework of the NLctM. In 
fact, one can start from the local expression of on the 
H-D line 



Ki-f/" COS {2VTrK$) + ... . (14) 



with K a constant and a the lattice spacing. In the same 
spirit as bosonization, the first uniform term comes di- 
rectly from Eq. (9). Since the Hamiltonian is quadratic in 
5*^ , the second staggered term would generate the first ir- 
relevant cosine operator (allowed by symmetries) that one 
could put in Eq. (13) at criticality. Now, independently 
of the scaling dimension of the exponential string in Eq. 
(2), the derivatives that appear as a consequence of sub- 
stituting Eq. (14) into the outer spins Sj and raise 
the overall scaling dimension by 1, and hence cannot be 
responsible for the numbers above. So, the leading order 
is obtained by retaining only the cosine term in the outer 
spin operators. Starting directly from the lattice formula- 
tion, it can be seen that the alternating term {—y cancels 



when combined with the exponential of the string, in such 
a way that the sum in Eq. (2) can be effectively taken only 
on the uniform part (defined as the average of the neigh- 
boring spins in the doubled lattice cell). So in the contin- 
uum limit we are simply left with J dxdx^ = ^{x) — ??(0). 
Summing up, we expect that the asymptotically dominant 
contribution to the string correlator comes from 



OUx) « (0|5t(0)<S(x)|0) , S{x) 



e'**^"' cos i — <P 
V9 



(15) 

and particularly from the terms exp(±i-^^), that have 

the sought decay exponent r]s,z = K/2. Note that the 
prefactor tt in the string of Eq. (2) is crucial for this. 

Along the same lines, if we move off-criticality the first 
dominant term in Eq. (14) becomes (— )^/" cos (V itK<1>) 
and, recalling that K = n/g, we see that Eq. (15) ac- 
quires a constant term due to the exponentials with op- 
posite arguments. We suspect that this is the mechanism 
according to which the string order parameters become 
nonzero in the Haldane phase. Rigorously speaking, one 
should compute the off-critical correlators using the sine- 
Gordon action and the outcome is expected to depend on 
the sign of /x. Here we can only notice that, semiclassically, 
when /z > the potential has a minimum in = while 
when /i < the minima lie at <P± = ±^Tr/4K. The dif- 
ference is that, in the first case cos {^/ttK(!>q) ^ while in 
the second one cos {^/'kK$±) = cos(±7r/2) = 0, which is 
the desired behavior. However, a systematization of these 
ideas is still underway. 



5 The XY2 phase 

We consider now the XY2 phase, which coincides with a 
region of small negative A and D < —2, getting narrower 
as D decreases. For large negative values of D we can 
resort to a perturbative study, outlined in the Appendix, 
which shows that the model is mapped onto an effective 
S=l/2 XXZ spin chain 



eff 



J 



L 

E 



(QX QX 
'->i 



sfsf 



+1 



ASfSi 



+1- 



(16) 



where J = l/\D\, A = 4A|I?|j- 1 and Si are spin-1/2 
operators (in particular = 2Sj). 

From the exact solution of the spin-1/2 XXZ model [5], 
we can therfore argue that our Hamiltonian, in the large 
\D\ limit, has a narrow critical region for— (2|£'|)^^<A< 
0, with elementary excitations carrying spin 1. For positive 
A the system is in an AFM phase, while for A < — (2|D|)~^ 
it has a ferromagnetic GS. We can conclude that, in the 
continuum limit, the system can be mapped effectively 
onto a Gaussian model, where the boson compactifica- 
tion radius and hence the critical exponents depend on 
A. Also, at every point of the X-D parameter space there 
is a horizontal direction along which a marginal opera- 
tor renormalizes the parameter K, and another direction 
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(A|D| = const., for D <C —1) along which the universahty 
class does not change. 

Inspired by the above theoretical results, we have per- 
formed a multi-target DMRG calculation at the point 
(A = —0.05, D = —5), which in our mapping corresponds 
to a S=l/2 XXZ model with A = 0, described by a free 
Dirac fermion with velocity u^*''^ = 1/\D\ = 0.2. Con- 
sistently with theoretical predictions, wc have found that 
only even spin sectors become gaplcss, while excitations 
with odd values of Sf^^ remain massive in the thermody- 
namic limit. Prom the scaling of the GS and of the first 
doublet, we have obtained Coo — —5.114607, = 
0.20752 ± 0.00006 and c^"*"") = 1.01991 ± 0.00003. Keep- 
ing in mind that here m has to be identified with twice the 
eigenvalue of S^^^ we extract K again from dio, now asso- 
ciated with the first level in S^^t = 2. We find K^^^'"'^ = 
0.9976 ± 0.0004, clearly compatible with the hypotesis of 
being at the FD point, as seen in Fig. 3 from the intersec- 
tion of the dimensions do±i and d±2o- 

As for the spin-spin correlation functions, it turns out 
that the transverse ones, which do not have a direct iden- 
tification in terms of spin-1/2 operators, have an exponen- 
tial decay. On the other hand, at the free Dirac point the 
longitudinal spin-1 correlation function is simply 4(iS'g>S'J) 
and sholud thearefore decay as (1 -I- , with ry^ = 

2 [12]. This alternating behavior is reproduced by DMRG 
calculations {L = 32,48,64,80,100 with M = 300) that 
yield the value r?^ = 2.06 ± 0.03. 



6 Conclusions 

In this paper we have studied the spin-1 AFM Hcisenberg 
chain with the inclusion of both Ising-like and single-ion 
anisotropy terms with the aim of investigating the critical 
properties of this model. In particular we have examined 
the massless phases that correspond to a c = 1 CFT. 

Our analysis starts from an analytical approach, aimed 
at the identification of the effective continuum field theory 
that describes the low-energy sector of the lattice model. 
We have found that all c = 1 phases are described by 
a free Gaussian model (with no orbifold construction), 
with continuously varying critical exponents. These re- 
sults have been then checked numerically, with the use 
of a multi-target DMRG algorithm, which allows for ac- 
curate calculation of many excited states together with 
string and ordinary correlation functions. The agreement 
between the two methods is remarkable. 

In particular, along the H-D transition line both the 
NLctM predictions and the DMRG data indicate that the 
Gaussian parameter K changes continuously from the value 
K = 2 (corresponding to the BKT point) for A = to 
K = 1/2 (corresponding to the SD point) for the tricriti- 
cal point where the H-D and the H-I transition lines meet. 
The latter result is at variance with some claims [7,10] 
according to which this point should be described by an 
SU(2)2 WZWN model. Moving away from the H-D fine, 
a gap opens up due to the relevance of the cosine oper- 
ator in Eq. (13). An analytical argument, supported by 
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Sb FD f[ BKT 

Fig. 3. Scaling dimensions of the three most relevant pri- 
mary operators, (m = ±l,n = 0), (m = 0,n = ±1) and 

(m = ±2, n = 0) as functions of K. We have indicated the 
three special values if = 1/2: SD, K ^ 1: FD; K = 2: BKT). 
The vertical dot-dashed line indicate the value of K = 0.9976 
corresponding to the point (A = —0.05, D = —5) (in the XY2 
phase) estimated through the first level dio (as indicated by 
the open circle). The dotted and dashed lines have the same 
meaning for the points on the H-D line {\ = 1,D = 0.99) 
and(A — 0.5, D = 0.65) respectively. In each case, the throe 
horizontal lines mark the values of the direct numerical esti- 
mates of the scaling dimensions. 



numerical estimates of the decay exponents, that relates 
the so-called twisting operator to the longitudinal string 
order correlation function is also skecthed. As far as the 
XYl phase is concerned, there are analytical and numer- 
ical evidences that the effective CFT is again a Gaussian 
model with K > 2. In this case the sine-Gordon operator 
is irrelevant. 

Finally, it is shown that, in the XY2 phase, the low- 
energy physics of the spin-1 model is equivalent to a XXZ 
spin-1/2 chain with an anisotropy parameter — 1 < zi < 1. 

We would like to thank L. Campos- Venuti, G. Morandi, S. 
Pasini and F. Ravanini for useful discussions. This work has 
been in part supported by the TMR network EUCLID (con- 
tract number: HPRN-CT-2002-00325). 



Appendix: Large Negative D 

For convenience, we consider the Hamiltonian (1) sub- 
tracting out a constant term DL. in order to study the 
D —00 limit, keeping A finite. At the zero order in 
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perturbation theory we simply have the Hamiltonian 



D\Y,{S^f + \D\L, 



(17) 



whose GS has degeneracy 2^, corresponding to all the spin 

configurations with = ±1, i.e. not containing zero's. 
This degeneration is lifted by the perturbation part 



(18) 



We denote with Pq and Pi the projectors onto the sub- 
spaces = {\si,S2, ■ ■ ■ ,sl) ■ Vi, \si) ^ |0)} and = 
{\si, S2, ■ ■ • , sl) ■ 3i, |sj) = |0)} respectively, so that we 
can write: 



[PiH^Po Pi{n° + y})Pi 



(19) 



We can now look for an effective Hamiltonian 'W^ de- 
scribing the low energy sector by projecting the resolvent 
operator Q{E) onto the subspace Z^ [11]: 

Pog{E)Po = Po[E - n]-'Po = [E- H'^^iE)]-^ . (20) 
We have 

= PoH^Po+PoH^Pi{Pi[E-{n°+n^)]Pir^PiH^Po 

(21) 

where we can consider the approximation in which E = 
and expand to the second order in H^. Since the off- 
diagonal part of Ti.^ connects Z° to Z^ creating a couple 
of zeroes, at the leading order we find 

H^^ff = PoH^Po - -J—PoH^p^n^Po 



L 

E 



2\D\ 



(22) 



1 



8|£»| 



This effective Hamiltonian is acting within Z° where we 
have only two local states per site. So, within this sub- 
space, we get an identification of the local spin-1 operators 
S" with spin-1/2 operators S" according to the following 
table. 



Spin 1 


Spin 1/2 




2S'| 




2SJ 
2S- 






2a+sA 


^1 





In terms of the new operators, we obtain an effective 
spin-1/2 XXZ model 



= ^ E (^f + Sphi + (4API + l)5|^+r) 

(23) 



N 
"MD\ 



where the sign in front of the x-y terms has been changed 
by means of the unitary transformation SJ {—ySJ, 
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